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runt. j Detailed current Edit ie 
: rative Part of 
- Standard. Single — 

S_R4& = Special AC ne 
MTH - Standard Single the sion Floating ACOSD 
MTHSACOSD. RG = Special ACOSD routine 
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MTHSACOS 3 Floating Point Arc-cosine routine 16-SEP-19 1:02:15 VAX/VMS M 4-0 P 1 
1-014 : eases Pa SE RESRTMACOS me 


-TITLE MTHSACOS caring Point Arc-cosine routine 
ACOS,ALOSD) 
IDENT /1-014/ File: MTHACOS.MAR Edit: JCW1014 


. 
FURR AAE EAA AERA EEA A AREER AEAAAAAAAAAAAEAAAAARAAARAAAAAATRETRETREAEE eee eee eee eee 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE eo THE 


NCLUSI : S_ SOF TW 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eORPORAT On NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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; ABSTRACT: 

; MTHSACOS is a function which returns the floating point arc-cosine in 
; radians of its single precision floating point argument. The call is 

3 standard call-by-reterence. 

; MTHSACOSD is a function which returns the floating point arc-cosine in 
; degrees of its single precision floating point argument. The call is 

; Standard call-by-reference. 


VERSION: 01 
HISTORY: 
; AUTHOR: 
Peter Yuo, 29-Jun-77: Version 01 


MODIFIED BY: 
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MTHRTL.SRCJMTHACOS .MAR; 1 (1) 


MTHSACOS ; Flo 
1-014 HISTORY ; Detailed Current Eds t History “SEP=- MTHRTL. 


-SBTTL HISTORY ; Detailed Current Edit History 


; ALGORITHMIC DIFFERENCES FROM FP=11/C ROUTINE: none 
3 Edit History for Version 01 of MTHSACOS 
Ons chang ed MOVF #0, RO into CLRL RO. 

SEERROR changed — SIGNAL. 


«es Change 
Chonged error handlin mechanism. Put error res 


Update conret ht notice. JBS 16-NOV-78 
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SRC tim THACOS .MAR; 1 


ult in RO (or RO:R1) 


before calling MTHSS$SIGNAL in order to allow user modify error result. 


-1979 


1-006 - 

1-009 = ¢ ~ roange, / ~~ I NVARG ,. «twa sage Bs, $9 -DEC-78 
1-008 - "to fhe PSECT d 5 TSnay=197 5 

a's 4 - De setae severe ls. ae Mayo 9 

1-010 - pee MTHSSQRT . SBL 27- sets? 

1-011 - ange JSB entry to MTHSACO RBG 28-Sept 

ee 3 ~ hdde degree entry p points. f RAN. 2 22-MAR-1981 

1-013 = Modified computation - x* for ixi>=1/2 

1-014 = Modified computation of - x*2 for 1/2<=iKXi< 


RNH g,Sept-8) 
=0.625 JCW 26-APR-83 
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1-014 Bns ; Declarative Part of Modul ae ets 91:98:17 MTHRTL.SRCJMTHACOS.MAR; 1 2 3) 


-SBTTL DECLARATIONS ; Declarative Part of Module 
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mm 

aon 
re 
ro 
zo 
> 
— 
a | 


; INCLUDE FILES: 


EXTERNAL SYMBOLS: 


BL 
TR MTHSSQRT_R3 
XTRN MTHSATAN-RS 


SATAND_R4 
RN MTHSK INVARGMAT 
TRN MTHSSSIGNAL | 
; EQUATED SYMBOLS: | 


SF_1.0 = *F1.0 i.1.0 

LF-PI_Ov_2 = *0007733a16 + 0940311 
LF_PI = *0007733a16 + *0940 1 
LF = *X000043B4 : 90 


90 
LF-180 *x00004434 : 180 
value 4 ; value.rf.r 


OF DB40C9 
OF DB4149 


44h 
00004434 
00000004 
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: 
00004080 0 
0 
0 
0 
Q 
0 
0 
0 


MACROS: none 
PSECT DECLARATIONS: 


-PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
; program section for math routines 


0000 


OWN STORAGE: none 


CONSTANTS: 
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gle Precis SEP=-1984 MTHRTL.SRCIMTHACOS.MAR; 1 
-SBTTL MTHSACOS = Standard Single Precision Floating ACOS 


3+ 

; FUNCTIONAL DESCRIPTION: 

ACOS = single precision floating point function 
ACOS(X) is computed as: 


C 
AC F 
t ) = ATAN(SQRT(1=X##2)/X). 

t ) = ATAN(SQRT((1=X)#(14#X))/X). 
X) = ATAN(SQRT(1=X*#2)/X) + PI. 
(X) = ATAN(SQRT((1=X)#(14X))/X) + PI. 


bs os hh 
ee i a ae 
—tt—-Oxx x 
"I 
A™™NA HHH 


CALLING SEQUENCE: 
ACOS.wf.v = MTHSACOS(x.rf.r) 
INPUT PARAMETERS: 


LONG = 4 : define coneeare multiplier 
x = 1 * LONG ; Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: floating arc-cosine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 


Signals: MTHS_INVARG if {Xi > 1 with reserved operand in RO (copied to 
the signal mechanism vector CHF$L_MCH RO/R1 by LIBSSIGNAL). 


Associated message is: ‘INVALID ARGUMENT’’. Result is reserved operand -0.0 


unless a user supplied (or any) error handler changes CHFS$L_MCH_RO/R1. 
NOTE: This procedure disables floating point underflow, enables integer 


; overflow. 
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401C 
6D 00000000°GF 9% 


~-ENTRY MTHSACOS, “M<IV, R2, R3, R4> 


3; standard call-by-reference entry 


; disable DV (and FU), enable IV 


MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


MOVAB G*MTHSSJACKET_HND, (FP) 


G 5 
ing Point Arc-cosine routine i as tb 91:9 :15 VAX/VMS Macro v04-00 Page 5 
S = Standard Single Precision Fl 6-SEP-1984 11:20:11 (CMTHRTL.SRCIMTHACOS.MAR;1 (4) 


set handler address to jacket 
handler 


MTHSACOS ; Floa 
1-014 MTHSAC 


case of an error in routine 

If an error, convert signal to user PC 
and resigna 

RO = 'X! = @value(AP) 

call special ACOS routine 

return with result in RO 


MOVF @value(AP), RO 
BSB MTHSACOS_R4 
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MTHSACOS Point Arc-cosine routine 16-SEP-1984 01: AX/VMS Macro V04-00 
1-014 Patties i", t 9] 98? MTHRTL. SR RCIMT THACOS.MAR; 1 


S_R4 = Special ACOS routine -SEP-1984 
! 188 -SBTTL MTHSACOS_R4 = Special ACOS routine 
! ! 23 ; Special ACOS - used by the standard routine, direct JSB call. 
"3 6 : CALLING SEQUENCE: 
1 187 ; save anything needed in RO:R4 
1 1 3 ; nate ; input in RO 
10 189; JSB MTHSACOS_R4 
} 139 : RSB ; return with result in RO 
1 192 ° 
8 1 138 MTHSACOS_R4: : poec tal ACOS routine 
1 194 MTHSACOS- RS ; Release 1 name 
54 50 $9 8 1 195 acy RO, R4 : save X in R4 
08 4 196 BNO TEST_FOR_1.0 : branch #7 ix! > 0 
01 138 § 
4 1) 3x20 
01 iz 
50 OFDB4OC9 BF 50 Ste 5 MOVF #LF_PI_OV_2, RO 3; RO = PI/2 
05 bis 02 RSB 3; return PI/2 if ix! = 0 
01D 05 ; 
1 ee rc 
01D 208 ° 
01D 09 TEST_FOR_1.0: 
50 8000 8F AA bp) 10 6icw PaaS 9 RO ; RO = IX! 
08 30 51 00 11 CMPF S*#5F_1.0 > compare {X! with 1.0 
0 18 00 iF BGEQ fr6_ 103 i : branch if !X! >= 1.0 
027 14 ; 
4 1g : 0 < ix! < 1.0 
057 «517° 
50 4020 8F 861 0027 18 CMPW #*x4020, RO ; Check for loss of significance in 
9c 14 002C 1 BGTR 1$ § computing 1 = x*2 and branch 
51 08 43 0 0 SUBF3 RO, #1, R1 ; RI =1- xX 
50 9 40 00 : 1 ADDF #1, RO ; RO=1+X 
50 1 44 00 § MULF Ri, RO ; RO = 1 = xX*2 
of 11 0038 BRB 23 3 Join normal flow 
50 0 44 A 4 18: MULF RO, ; RO = Xee2 
58 $8 50 43 D 5 SUBF RO! § F_1.0, RO 3 RO = 1.0 = Xe02 
OOO0000'EF 16 41 § 23: JSB MTHSS Tas” ; RO = Sant Cloasss) 
50 5 46 47 IVF Ran R RON ; RO = SQRT(1-Xe*" 4 
4 oD 4A 3 PUSHL R 3 save sin 2 $97 test 
00000000" EF 16 4 JSB MTHSATAN_R4 3; RO = AT Me SaRT Cloneee /x) 
54 8E€ D0 005 0 MOVL (SP)+, R& 3 restore sign of xX 
046 11 09 1 BRB TEST_SIGN > branch to TEST_SIGN 
bs : : 
ef : : 1 s< ix! 
? a” 
3f $ GEQ_T0_1.0: 
OE 14 57 8 BGTR ERROR 3; branch to ERROR if ‘xX! > 1.0 


=o 
. 
MTHSACOS ; Floating Point Arc-cosine routine 16-SEP-1984 01:02:15 VAX/VMS Macro v04-00 Page 7 | 
1-014 ATHSAcOS. RG = Special ACOS routine aoe Pi 945981} EMTHRTL. SRENTHACOS..MAR: 1 5) | 
39 
4 
41; 
9 62 3 Iki s 1.0 
9 43; 
23 44 | 
50 04 5 45 CLRF RO ; RO = 0 
58 é$ 
058 4 
058 ry 3 
bee rts ; Test the sign of X in order to decide if add FI to the result 
beB 584 
058 26 TEST_SIGN: 
54 53 0058 5 TSTF R4 ; test the sign of xX 
07 18 0050 54 BGEQ 10$ ; branch if X > 0 
50 OFDB4149 8F 40 OO0SF 55 ADDF WLF_PI, RO ; add PI to RO if x <0 
05 966 2$ 10$: RSB 3 return with result in RO 
06 5 
067 58 ; 
bee) 59; 1 < IX!, error 
067 60 ; 
0067 61 
6— 0D 0067 6¢ ERROR: PUSHL (SP) 3: return PC from JSB routine 
7E OO°8F 9A 0069 6 MOVZBL #MTHSK_INVARGMAT, -(SP) ; condition value 
50 01 OF 78 006D 64 ASHL #15, #T, RO ; RO = result = reserved operand -0.0 
865) 65 3 pres to signal mechanism vector 
071 66 ; CCHFSL_MCH_RO/R1) so error handler 
a 67 3 can modify the result. 
00000000'GF 02 FB 0071 68 CALLS #2, G*MTHS$SIGNAL : signal error and use real user's PC 
0078 69 ; independent of CALL vs JSB 
05 Bne8 4 RSB ; return - RO restored from CHF$L_MCH_RO/R1 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: floating arc-cosine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signals: MTHS_INVARG if iX! > 1 with reserved operand in RO (co180ed to 
the signal mechanism vector CHFSL_MCH_RO/R1 by LIBSSIGNAL). 


Associated message is: "INVALID ARGUMENT'’. Result is reserved operand -0.0 
unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1. 
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MTHSACOS Float Point Arc-cosine routine 16-SEP-1984 01:02:15 VAX/VMS Macro V04-00 Pa 8 
1-014 MTHSACO . = Standard Single Precision F moet 9139 711 EMTHRIL. SREIMTHACOS MAR: 1 ” (6) | 
4 -SBTTL MTHSACOSD = Standard Single Precision Floating ACOSD 
7 
76 344 
4 ; FUNCTIONAL DESCRIPTION: 
7 ; ACOSD = single precision floating point function 
: ACOSD(X) is computed as: 
; If X = 0, then ACOSD(X) = 90. 
3 If X = 1, then ACOSD(X) = 0. 
; If xX = <1, then ACOSD(X) = 180. 
: If O<x< V2. then ACOSD(X) = ATAND(SQRT(1-X*#2)/X). 
; If 1/2 <x< then ACOSD(X) = ATAND(SQRT((1-X)#(14x))/xX). 
; If -1/2 < x < 6, then ACOSD(X) = ATAND(SQRT(1-X##2)/XxX) + 180. 
89 ; If -1 < X < =1/8, then ACOSD(X) = ATAND(SQRT((1-X)#(14x))/X) + 180. 
90 ; If 1< iXt, error. 
: CALLING SEQUENCE: 
; ACOSD.wf.v = MTHSACOSD(x.rf.r) 
0 : INPUT PARAMETERS: 
00000004 LONG = 4 : define longword multiplier 
00000004 p x = 1 * LONG ; Contents of x is the rement 
0 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


Se Se Ge Ge Ge Ge Ge Ge Fe Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge 


-ENTRY MTHSACOSD, “M<IV, R2, R3, R4> 

3 standard call-by-reference entry 

: disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


NAME WIN 9 OD NA UE WIN 9 OD NAME WIN SO OD NA NE WIN = O OD NAME WN OOONOUSW 


nn > > > > > > > > > > > > > > > > > >> b> > > > > > > > > > >> > > > > > > > > > > > > > > > > > > > > i ile 
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401C 
6D 00000000°GF 9E 


MOVAB G*MTHSSJACKET_HND, (FP) 


ngle ariiten F 


ewnr—O Oe 


715  VAX/VMS Macro v04-00 Page 
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; set handler address to jacket 
3; handler 


case of an error in routine 

If an error, convert signal to user PC 
and resigna 

RO = iX! = *avalue(A P) 

call special ACOSD routine 

return with result in RO 


o V04-00 Page 
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R4 = Special ACOSD routine -SEP-1984 MTHRIL. 
-SBTTL MTHSACOSD_R4 = Special ACOSD routine 
; Special ACOSD - used by the standard routine, direct JSB call. 


+ CALLING SEQUENCE: 
save anything needed in RO:R4 


eee R 
JSB MTHSACOSD_R4 
RSB 


MTHSACOS 
1-014 


; Flo 
MTHSA 


3; input in RO 
3 return with result in RO 


MTHSACOSD_R4:: ; special ACOSD routine 
54 50 39 MOVE RO, R4 i; save X in R4 
08 BNEQ  D_fEST_FOR_1.0 : branch if {x! > 0 
> Kk =0 
50 00004384 8F 50 MOVF #LF_90, RO ; RO = 90 
05 RSB 3 return 90 if ix! = 0 


* ¥*x8000, RO ; RO = IX! 
CMPF =—s«RO,-S*#SF_1.0 compare !X! with 


50 8000 
BGEQ 0_GEQ_T0_T.6 branch if |X! >= 


co 
Wwwoo 
oon 
—w 
o> 
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oo 


:0< IX! < 1.0 
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50 4020 8F 81 CMPW #*x4020, RO ; Check for possible loss of 
OC 14 GT 1$ § significance with 0.625 
51 08 50 43 SUBF3 RO, #1, R1 ; RI =Te- xX 
50 o¢ 40 ADDF #1, RO ; RO=1 + X 
50 1 44 A MULF R1, RO ; RO = 1 = X*2 
o? 11 ONB BRB 23 3; Join main flow 
50 0 44 008 1$ MULF RO, RO ; RO = X*#2 
50 08 50 43 008 SUBF RO, S*#SF_1.0, RO 3 RO = 1.0 = X#*2 
90009300 "5! 16 23: JSB MTHSSQRT_RS 3; RO = hb Bt tf 
5 4 46 DIVF R4, RO 3 RO = SQRT(1-X*#2)/xX | 
54 DD PUSHL R 3 save sign of X for sign test 
00000000 ‘EF 16 JSB MTHSATAND_R4 ; RO = ATAND(SQRT(1-X**2)/X) 
54 aS dO MOVL (SP)+, R4 3; restore sign of X 
07. «11 BRB D_TEST_SIGN : branch to TEST_SIGN 


: branch to ERROR if ‘xX! > 1.0 


OOOO 090909 69 G9 G9. G9 GD G0 i NDS DS DDD DDO TTI BB ES EE EWI 
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MTHSACOS ; Floating Point are 19-36 0= 1388 91:98:18 yaa lyase Macro V04-00 Page 11 
1-014 ATH ACOSD_R4 = Special ACOSD routine -SEP-1984 11:20:11 CMTHRTL.SRCIMTHACOS.MAR; 1 (7) 
D 93 
D 8 3 
D 5 3; ixt = 1.0 
D8 397 
50 D4 4 38 10$: CLRF RO ; RO = 0 
ye 
O07 106 : Test the sign of X in order to decide if add 180 to the result 
Ob? 404 : 
D? «405 D_TEST_SIGN: 
54 53 O00D7 406 TSTF R4 ; test the sign of X 
tH 18 Ba08 407 BGEQ 10$ ; branch if xX > 
50 00004434 8F 40 00D 408 ADDF #LF_180, RO ; add 180 to RO if x < 0 
05 Boe +4 10$: RSB 3; return with result in RO 
bor 411 - END 


i kj Sat 9 91:06:13 tT Macro V04-00 
6-SEP=1984 11:20:11 CMTHRTL.SRCIMTHACOS.MAR; 1 


N 
MTHSACOS ; Floating Point Arc-cosine routine 
Symbol table 
GEQ_1T0.1.0 000000 R 1 
D-TEST_FOR 1.0 09 009 R 1 
D-TEST_SIGA D7 R 1 
ERROR 0 Bh R 4 
GEQ_10_1.0 0000057 R 1 
LF_T80 = 00004434 
LF_90 = 00004384 
LF_PI = OFDB4149 
LF _PI_OV_2 = OFDB40C 
G = 00000004 
alba gy ody HND *eeeetes x 1 
MTHSSSIGNA eeneeeee 0 
MTHSACOS 00000000 RG 1 
ACOS 0000079 RG 01 
MTHSACOSD_R4 0000089 RG 01 
MTHSACOS_R4 B80 0010 RG 01 
MTHSACOS RS 00010 RG 01 
MTHSATAND_R4 tereeree§ =X 00 
MTHSATAN_R4 etetetee x 00 
MTHSK_INJARGMAT teeeeene x 00 
MTHSSGRT_R3 teeerere =X 00 
SF_1.0 = 00004080 
TEST_FOR 1.0 0000001D R 01 
TEST SIGN 00000058 R 01 
ALUE = 00900004 
ee ee + 
! ; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
ABS 00000000 < 9.3 89 0.) NOPIC USR CON ABS LCL N 
MTHSCODE OOOOO0ES ¢ 227.) te PIC USR CON REL LCL 
neeetebneeemeceneunsaadand 
H ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 30 0: 8:00.97 0:00:01.11 
Command processing 115 0:00:00.59 0:00:03.37 
Pass 1 94 :00:01.1 8:88: 9 
syabol iable sort 5 :00: 0:3 :00: 4 
wt table output 0:00: 03 0:00:00 08 
Psect synopsis output 0: 8: 0 0:00:00.1 
Cross-reference output 0 0: 2 -00 Ba RR ay 00 
Assembler run totals 329 00:00:02.80 :00:11.48 


Limit was 700 pages. 


The working is 
f pages) of virtual memory were used to buffer the intermediate 


6148 bytes 
There were 
471 source (Snbs were read in Pass 1, producing 


object records in Pass 
1 page of virtual memory was used to “define 1 macro. 


ode. 
0 pages of symbol table space allocated to hold 26 non-local and ; local symbols. 


Page 


OSHR NOEXE NORD NOWRT NOVEC BYTE 
SHR EXE RD NOWRT NOVEC LONG 
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$e omen enema crew n some seneeee =} 


! Macro Library statistics ! 


¢eeeeeeeoceecee waemamacao sac a } 


Macro Library name Macros defined 
_$255SDUAZ8: (SYSLIBISTARLET.MLB;2 ta at We 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHACOS/OBJ=OBJ$:MTHACOS MSRC$:MTHJACKET/UPDATE=(ENH$:MTHJACKET) #MSRCS: 


0257 vaxvums va ENT CORPORATION 


V4.0 3 . ND PROPRIETARY | 


